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We have developed an analytical model based on the perturbation theory in order to study the 
optical propagation of two successive intense solitons in hollow-core photonic crystal fibers filled 
with Raman-active gases. Based on the time delay between the two solitons, we have found that 
the trailing soliton dynamics can experience unusual nonlinear phenomena such as spectral and 
temporal soliton oscillations and transport towards the leading soliton. The overall dynamics can 
lead to a spatiotemporal modulation of the refractive index with a uniform temporal period and a 
uniform or chirped spatial period. 


Raman scattering involves the transfer of a small frac¬ 
tion of energy of an optical field to another field at longer 
wavelengths via inelastic scattering with optical phonons 
supplied by the medium. The Raman self-frequency red- 
shift, first observed in conventional optical fibers [H-Q; 
can continuously downshift the central frequency of a 
single pulse during propagation. With the invention of 
solid-silica-core photonic crystal fibers (PCFs) ii , this 
process has been demonstrated and extensively exploited 
in attaining a broad supercontinuum Q. The main in¬ 
gredients are the high confinement of light and tunability 
dispersion provided by these fibers. Gas-filled hollow- 
core (HC) PCFs with Kagome lattice have become a 
strong competitor for solid-silica-core PCFs 0,0. Be¬ 
sides the latter benefits, a wide range of gases with differ¬ 
ent properties offers several opportunities to demonstrate 
new phenomena and applications using optical fibers 0 - 
0 . 

Raman scattering processes in gases are characterized 
by having very long molecular coherence relaxation (de¬ 
phasing) time (~ps), at least three orders of magnitude 
higher than silica glass. Within this time, the medium 
possesses a high non-instantaneous nonlinear response. 
A Raman-induced continuous downshift of the frequency 
of an ultrashort pulse propagating in a gas has been ob¬ 
served [IJ. Moreover, a sinusoidal temporal modulation 
of the medium refractive index lagging the pulse has been 
detected by launchii^a delayed weak probe within the 
dephasing time (2ll - [^ . Very recently, we have treated 
this modulation as a periodic temporal crystal [2^ and 
we have proposed to use hydrogen-filled HC-PCFs to ob¬ 
serve the analogue condensed matter physics phenomena, 
such as the Wannier-Stark ladde r , the Bloch oscilla¬ 
tions and Zener tunneling |26j|. In this Letter, we 
investigate the case when the delayed pulse is another 
strong fundamental soliton rather than a weak probe 
pulse as in [2^. The delayed soliton can be treated as 
a particle in a moving periodic potential induced by the 
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leading soliton. Phenomena such as spectral and tempo¬ 
ral soliton oscillations as well as upgrading the temporal 
crystal to a spatiotemporal one are predicted here. 

Pulse propagation in Raman-active gases is governed 
by the interplay between the Maxwell equations and the 
Bloch equations of an effective two-level system describ¬ 
ing a single Raman transition. We have previously shown 
that this set of equations can be simplified to a single gen¬ 
eralized nonlinear Schrodinger equation for femtosecond 
pulses with an energy of a few micro-Joules [^ . 

= 0 , ( 1 ) 

where i? (r) = k siii[S {t — T')]\'ijj {t')^ dr' is the 
nonlinear Raman convolution, ip is the pulse complex en¬ 
velope, ^ is the longitudinal coordinate along the fiber, r 
is the time variable in a reference frame moving with the 
pulse group velocity, S = 27 r/A is the Raman frequency 
shift, A is the corresponding oscillation period, k is the 
relative strength between the Raman and Kerr nonlinear¬ 
ities, and the equation is in normalized units. We have 
assumed propagation in the deep anomalous dispersion 
regime to determine the pure effect of the temporal non¬ 
locality of the Raman nonlinearity on successive intense 
solitons, which can occur during the soliton fission pro¬ 
cess in supercontinuum generation in gases [l^ . 

We are interested in investigating the propagation of 
two successive strong ultrashort pulses, ip = V'l + ' 02 - 
The two pulses are assumed to have the same frequency, 
hence, they will propagate with the same group velocity, 
and experience the same dispersion. Also, they are not 
temporally overlapped, but separated by a delay less than 
the Raman dephasing (relaxation) times. In this case, 
all the nonlinear cross terms are canceled out, with the 
exception of the Raman polarization effect of the leading 
pulse on the trailing pulse. Under these conditions Eq. 
(HD can be split into two equations for ipi and 02, 

+ \ipj + Rj {t) ipj = 0, (2) 

with j = 1,2. For weak Raman nonlineari¬ 

ties, the solutions of these equations are two 
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perturbed fundamental solitons, ipj (^, r) = 
Vj sech \Vj (r — fj (^))] exp [—i^j (^) r] where Vj, 
and fj are the soliton amplitude, central frequency, 
and temporal peak, respectively. When the soliton 
durations <C 1/(5, their Raman response functions can 
be approximated by using a Taylor expansion [ 2 ^, 

Rj ('r) ~ X! ("T - ^i)]} ! 

i 

( 3 ) 

with I = 1,2. Hence, the leading soliton will induce a 
retarded sinusoidal Raman polarization that will impact 
the dynamics of the trailing soliton. 

To understand how Raman nonlinearities can affect the 
pulse dynamics, we have adopted the variational pertur¬ 
bation method [l^. Using this method, we have con¬ 
sidered each Raman contribution Rj as a perturbation 
for the associated solution (a fundamental soliton j) of 
its equation. A set of coupled governing equations that 
determine the evolution of each soliton parameters can 
then be derived, 

Hi = —91 

U =giC^/2, 

H 2 = - 92 ^- 92 ^ [ cos [5 (f 2 - fi)] d^, (4) 

^ k2 Jo 

T2 = - f H 2 
Jo 

where gj = i«;7r(5^csch (7r(5/2Vj). Hence, the first (lead¬ 
ing) soliton will always linearly redshift in the frequency 
domain with rate gi , and decelerate in the time domain. 
Whereas the dynamics of the second (trailing) soliton 
depends on two different components: (i) its own (self) 
component that will lead to a linear redshift similar to 
the leading soliton, with rate g 2 ', (ii) a cross component 
representing the effect of the first soliton on the second 
soliton. The latter component is proportional to the ra¬ 
tio between their amplitudes and the cosine of the time 
difference between them. Since the cosine term varies be¬ 
tween positive and negative values, the dynamics of the 
second soliton can switch back and forth between redshift 
and blueshift in the frequency domain or deceleration and 
acceleration in the time domain. This analytical model 
shows a very good agreement with the numerical model 
as depicted in Fig. [U provided that the assumption of 
the soliton durations ^ 1/(5 is satisfied. Deviations from 
this condition results in an overestimated value of the 
frequency-shift rate, g, of an about a factor of 2. In Fig. 
[U t = Tto, z = ^zq are the physical delay and propaga¬ 
tion distance, respectively, zq = tQ/|/52|, to is the pulse 
duration, /32 is the second-order dispersion coefficient. 
Also, this method can be extended easily to study the 
interactions of more than two solitons. 

Figured] shows four special cases of the temporal evo¬ 
lution of the second soliton superimposed on the total 
induced-Raman polarization in a reference frame moving 
with the leading soliton t = t — g\^ V 2 is chosen 
smaller than Vi so that the cross component in Eq. Sjis 



FIG. 1. (Color online). (a,b) Spectral and temporal evolution 
of two fundamental solitons launched in a HC-PCF filled by 
a Raman-active gas, respectively, with k = 0.16, and S = 0.5 
(equivalent to a Raman-mode oscillation with period 185 fs, 
such as in deuterium [^1. All subsequent simulations in 
this Letter use these parameters. The leading and trailing 
solitons have normalized amplitudes Ui = 2.5, V 2 = 1.25, and 
full width half maximum (FWHM) 8 fs, 16 fs, respectively. 
The black dashed-dotted lines are the analytical predictions 
using Eq. Contour plots are given in a logarithmic scale 
and truncated at -40 dB. 
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FIG. 2. (Color online). Raman polarization induced by two 
fundamental solitons propagating in the gas-filled HC-PCF, 
used in Fig. [T1 in a reference frame moving with the leading 
soliton. White solid lines represent the temporal evolution of 
the two solitons using the analytical prediction [Eq. ([4|]. The 
second soliton is launched at (a) f 2 (0) = 6.75A, (b) f 2 (0) = 
6.5A, (c) f 2 (0) = 5 . 25 A, and (d) f 2 (0) = 5A. 


comparable to the self component. The leading soliton 
induces a lagging sinusoidal temporal modulation of the 
medium refractive index due to the Raman effect, result¬ 
ing in a temporal crystal with period A. Due to soliton 
acceleration or spectral redshift, a constant force —gi 
acts on the crystal in the positive-delay (right) direction 
[^ . Hence in the accelerated moving frame, there is an 
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effective force acting on the trailing soliton on the other 
direction, similar to a particle in a moving periodic poten¬ 
tial. The uniformity of the crystal is modified along the 
direction of propagation, resulting in a spatiotemporal 
modulation of the refractive index, i.e. a spatiotemporal 
crystal . As we are operating in the anomalous disper¬ 
sion regime, the positive (negative) variation of the re¬ 
fractive index represents a potential well (barrier). Based 
on the initial time delay between the two solitons Afv, 
the dynamics of the second soliton is different, and the 
shape of the crystal is altered. Looking at Fig. [51 launch¬ 
ing the second soliton at (a) the top of a barrier or (b) 
at the right edge of a well induced by the first soliton, 
the second soliton will be able to overcome the barriers 
during propagation, and transported across the poten¬ 
tial by the applied force to the left direction. The output 
spatiotemporal crystals are chirped along the direction 
of propagation in these cases. Interestingly, the second 
soliton in (b) experiences a net maximum self-frequency 
blueshift of 9.12 THz = 51.4 nm after 10 cm of propaga¬ 
tion, before it is redshifted. Launching the second soliton 
at (c) the minimum or (d) a left edge of a well, the second 
soliton will not be able to overcome the barriers in these 
cases, so it is trapped inside the well and will oscillate 
indefinitely. The amplitude of oscillation in (c) is very 
small, since the initial velocity of the soliton in this po¬ 
tential is zero. The soliton will oscillate in an asymmetric 
manner due to the modified potential beyond the second 
soliton peak as well as the acting force that is opposite 
to the initial velocity as in (d). The resulting spatiotem¬ 
poral crystals have uniform periods along the direction 
of propagation in both cases. 

The temporal evolution of two successive ultrashort 
Gaussian pulses rather than fundamental solitons are 
portrayed in Fig. |3] for different time delays. The two 
pulses have the same central frequencies and amplitudes, 
and the delay is within the relaxation coherence time 
of the Raman-active gas. The two pulses will experience 
pulse compression and soliton fission processes. The lead¬ 
ing ‘first’ pulse excites the Raman polarization ‘potential’ 
that will affect the trailing ‘second’ pulse. The dynam¬ 
ics of the leading pulse is independent of the time delay 
and will encounter a self-induced Raman redshift (decel¬ 
eration). The trailing pulse dynamics is influenced by 
its Raman-induced effect as well as the cross-accelerated 
Raman polarization effect of the leading pulse. In Fig. [31 
panel (a) shows the case when the self and cross compo¬ 
nents are working together, resulting in a strong delay in 
comparison to the first pulse. The situation when the self 
and cross components acts against each other is figured 
in panel (b), where initially the second pulse is nearly 
halted since the cross and self components cancel each 
other. When the second pulse is launched at one min¬ 
ima of the potential induced by the first pulse, the pulse 
is well-confined during propagation, even after the pulse 
fission the generated solitons are still traveling together, 
see panel (c). Launching at one potential-maxima, the 
dynamics of a tree-like behavior is obtained as shown in 



FIG. 3. (Color online). Temporal evolution of two successive 
identical Gaussian pulses with profile A exp (—t^/ 2) in the 
gas-filled HC PCF, with N = 7, and FWHM = 25 fs. The first 
pulse is launched at ri (0) = 0. The second pulse is launched 
at: (a) f 2 (0) = 2A, (b) f 2 (0) = 2.5A, (c) f 2 (0) = 2.25A, and 
(d) r2 (0) = 2.75A. 


panel (d), where each soliton propagates in a different 
direction. In all the cases, the dynamics of each soliton 
depends on where exactly this soliton is emitted inside 
the total accelerated periodic potential induced by other 
preceding solitons. 

In conclusion, the propagation of two successive ultra- 
short strong pulses in HC-PCFs filled by Raman-active 
gases have been investigated. The two pulses are not 
temporally overlapped and separated by a delay smaller 
than the Raman polarization dephasing time (~ 100 ps). 
We have adopted perturbation theory to understand the 
effect of the Raman nonlinear polarization on the dynam¬ 
ics of the two pulses, assuming that they are fundamental 
solitons with durations much less than the Raman tran¬ 
sition period. We have found that the leading soliton will 
always redshift, while the trailing soliton dynamics can 
oscillate between redshift and blueshift processes based 
on the time delay between the two solitons. The sec¬ 
ond soliton can experience a boosted blueshift without 
experiencing losses, in contrast to the ionization-induced 
blueshift technique MM- The leading soliton induces 
an accelerated temporal crystal due to the long-living 
Raman polarization. This crystal is upgraded to a spa¬ 
tiotemporal one with a uniform or chirped spatial period 
beyond the trailing soliton. Our simulations have also 
been extended to include the case of two energetic Gaus¬ 
sian pulses. We have found that a good control of the 
trailing pulse dynamics can be achieved via changing the 
time delay between the two pulses. These theoretical re¬ 
sults will pave the way to manipulate and control pulse 
dynamics in PCFs for novel optical applications. 
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